(N 



CD 

(N 
> 



X 



INHA-NTG-Ol/2012 

Parity-violating vrA^A^ coupling constant in the chiral quark-soliton model 

Hee-Jung Lee/'B Chang Ho Hyun,2^3'|t] and Hyun-Chul Kim4.5.3.| 

' Departament of Physics Education, Chungbuk National University, Cheongju 361-763, Republic of Korea 

^Department of Physics Education, Daegu University, Gyeongsan 712-714^, Republic of Korea 

'^School of Physics, Korea Institute for Advanced Study, Seoul 130-722, Republic of Korea 

'^Department of Physics, Inha University, Incheon ^02-751, Republic of Korea 

^Department of Physics, University of Connecticut, Storrs, CT 06269, U.S.A. 

(Dated: June 2012) 

We investigate the parity- violating ttNN Yukawa couphng constant h\f^f^ within the framework of 
the SU(2) chiral quark-soliton model, based on the AS' = effective weak Lagrangian derived within 
the same framework. We find that the parity-violating tvNN coupling constant is about 1 x 10~ 
at the scale of 1 GeV. The results of ftJriViv turn out to be sensitive to the Wilson coefficient. We 
discuss how the gluonic renormalization supresses the parity- violating ttNN coupling constant. 
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I. INTRODUCTION 

Oj' The parity-violating (PV) hadronic processes in low-energy regions have been one of the most fundamental issues 

I in nuclear and hadronic physics for long time (see a recent review [l| for some historical and phenomenological 

j^ background). However, the weak interactions of hadrons are yet poorly understood because of the strong interaction, 

compared to lepton-lepton or lepton-hadron weak processes. For example, the long-standing puzzle of the A/ = 1/2 

rule in strangeness-changing weak interactions indicates that the effect of the strong interaction in weak processes 

raises a non-trivial problem [2-13|- It is even more difficult to study parity- violating nuclear processes because of 

experimental feasibility and theoretical complication caused by the nonperturbative strong interaction of quarks and 

gluons. The standard model (SM) asserts that charged weak boson exchange induces flavor-changing weak interactions 

^j-^ ' whereas the neutral current conserves the flavor. The basic ingredient to describe low-energy hadronic weak processes 

h — is the quark current-current interaction with W and Z bosons. However, in order to describe low-energy phenomena 

■^ , below 1 GeV, one has to scale down this interaction from the mass scale of the W and Z. In the course of this scaling, 

the quark-gluon interactions are encoded in the Wilson coefficients by the renormalization group equation (5|-l8| , which, 

(^ ', however, explains only a perturbative part of the strong interaction. 

p\| ' Desplanques, Donoghue, and Holstein (DDH) ^ suggested that hadronic and nuclear PV processes can be described 

by one-boson exchange such as n-, p-, and w-exchanges [8l-[l0l| a la the strong nucleon-nucleon (NN) potential. The 
main factors of the PV NN potential are the seven weak meson- iViV coupling constants, i.e. h\jqj^, ^^pNN^ ^]>nn^ 
^'pNN^ ^uiNN^ ^IjNNj ^^'^ ^'pNN^ where superscripts denote the isospin difference A/. Among these coupling constants, 
it is of utmost importance to understand the PV ttNN coupling constant, because it governs the long-range part of 
^ I the PV NN interaction, so that it plays the most signicant role in explaining the PV nuclear processes. The PV ttNN 
■ ■ coupling constant can in principle be extracted from various PV reactions np — )■ d'j |ll| - [l3j . and ^^F* — >-^^ F [ij- 
Il6| but is fraught with large undertainties. We refer to a recent review '17| for the present status of hadronic PV 
experiments. It has been also calculated in various theoretical frameworks: the SU(6)^y quark model [allMl with the 
effective weak Hamiltonian, the Skyrme model [19l - [2]| . and QCD sum rules ^2^, and so on. Even though a great 
amount of efforts was made on understanding /i^jvAf experimentally as well as theoretically, its quantitative value is 
still elusive. 

In the present work, we investigate the PV ttNN coupling constant, h\jqj^, within the framework of the SU(2) 
chiral quark-soliton model (xQSM) which is an effective chiral model for QCD in the low-energy region with con- 
stituent quarks and the pseudoscalar mesons as the relevant degrees of freedom. The model respects the spontaneous 
breakdown of chiral symmetry and describes baryons fully relativistically. Moreover, it is deeply related to the QCD 
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vacuum based on instantons [23| and contains only a few free parameters. These parameters can mostly be fixed 
to the meson masses and meson decay constants in the mesonic sector. The only remaining free parameter is the 
constituent quark mass or dynamical quark mass that is also fixed by reproducing the electric properties of the proton. 
The xQSM was successful in describing lowest- lying baryon properties [2J|. Furthermore, the renormalization scale 
for the xQSM is naturally given by the cut-off parameter for the regularization which is about 0.36 GeV^. Note that 
it is implicitly related to the inverse of the size of instantons {p « 0.35 fm) [25|, [26|. This renormalization scale is very 
important in general, because the essential feature of the PV hadronic interactions comes from the effective weak 
Hamiltonian that has a specific scale dependence, as mentioned previously. Thus, the matching of this scale consists 
of an essential part in investigatng any nonleptonic decays and PV hadronic processes. 

While the xQSM provides a plausible framework to study the PV nNN coupling constant, there are at least 
two theoretical difficulties. Firstly, the effective weak Hamiltonian has two-body operators and one has to treat the 
four-point correlation functions in order to compute the PV ttNN coupling constant. Secondly, since the momentum- 
dependent dynamical quark mass is known to play a significant role in describing K ^ tttt nonleptonic decays [27|, |28| , 
one can expect that it would also contribute to h\j^j^ substantially. This is in particular important, because a certain 
amount of non-perturbative effects is refiected in the momentum-dependent quark mass, which arises from the zero 
mode of instantons. However, it is very difficult to handle these problems in the self-consistent xQSM. In order to 
circumvent all technical difficulties in the self-consistent approach, we will use the gradient expansion to calculate 
^ttNN^ taking the limit of a large soliton size, so that valence quarks in a nucleon plunge into the Dirac sea and the 
soliton emerges as a topological one J24| . which is quite similar to a skyrmion. Equivalently, we can start directly 
from the AS' = effective weak chiral Lagrangian derived in Ref. [29| and quantize the chiral soliton collectively. 
Then, we introduce a physical pion through quantum fluctuations around the soliton field. This procedure will lead 
to the results for /i^^vAf without fitting any parameter. In the present work, we will restrict ourselves the SU(2) case 
for simplicity and will concentrate on how the low-energy constants (LECs) found in Ref. [20| feature the PV ttNN 
coupling constant. 

This paper is organized as follows: In Section [Hi we describe briefly a general formalism for the derivation of the 
PV weak ttNN coupling constant. In Section Hill we present the numerical results for h-^j^j^ and discuss the role of 
the LECs of the AS = effective weak chiral Lagrangian. The last section is devoted to the summary and outlook of 
this work. 

II. GENERAL FORMALISM 

In this Section, we will show how to incorporate the AS = effective weak Hamiltonian into the effective chiral 
action. We employ the AS — effective weak Hamiltonian derived in Ref. Q. The Hamiltonian reads 
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where the operator 0{Mi,Ni) is defined as a two-body operator 0{Mi,Ni) = —ip^'ffj^jc^Miipip^j^Niip in Euclidean 
space, and tA denotes the generator of the color SU{3) group, normalized as tri^is = 25 ab- The definitions of the 
matrices Ai and Bi, and the coefficients a, /3, 7 and p can be found in [8]. These coefficients are the functions of the 
scale-dependent Wilson coefficient K{p) defined as 






Kip)^ll + '-^b\n^\, (2) 



where g{p^) denotes the strong running coupling constant, p stands for the renormalization point that specifies the 
energy scale, 6=11 — 27V//3, and M\y is the mass of the W boson. The coefficient K encodes the effect of the strong 
interaction from perturbative gluon exchanges. 

The four-quark operators are expressed generically by 

Q\x) = -i^Hx)T\i;ix)i;HxT2^i^) , (3) 

where z(= 1, • • • , 12) labels each four-quark operator in the effective weak Hamiltonian and ^\(2) consist of the Dirac 



gamma and flavor matrices. Thus, the effective weak Hamiltonian can be rewritten as follows: 



12 



<^=" = ^C.Q'(x), (4) 



j=i 



where Ct denotes a, /3, 7 and p according to Eq. ([TJ. 

In order to derive h^j^j^ in the xQSM, we have to solve the following matrix element: 



12 12 



The nucleon state is defined in terms of the loffe-type current in Euchdean space {xq — —ix4): 

\N{pi)) = lim eP^y^Af*ipi) ( d'ye'P^y f^{y)\Q), {N{p2)\ = Jiin_ e-P°"W(p2) f d^xe-'P^-^{0\JNix). (6) 



2:4— f 4- 00 

The nucleon current Jj^ (Jn) pla-ys a role of creating (annihilating) nucleons. The J\f* [J\f) represensts the normalizing 
factor depending on the initial (final) momentum. The JJ^ (Jn) consists of Nc quarks: 

Jn{x) = ^e^^^^--""^r^T47KJj3F«)^«ic,(a:)---^«„,c„Ja:), (7) 

where si • • • sn^ and ci • • • cn^ denote respectively spin-isospin and color indices. The r}^,^ y)( u y ) ^^^ matrices 
with the quantum numbers {TT3Y){JJ3Yi^). For the nucleon, T = 1/2, F = 1 and J = 1/2. The right hypercharge 
will be constrained by the baryon number. The creation baryon current is written as 

jUv) - ]^^""-^"=r--^-;;^^) (-^^'74),„^,„^ (X) . . • (-.^^74),,,, (x). (8) 

The partial conservation of the axial- vector current (PC AC) being considered, the matrix elements in Eq.(l5|) can be 
related to the following four-point correlation function 

12 
lim J2c,{0\T[JN{x)Q\z)d^Al{OjUy)m = lim ^, (9) 

XQ^- + cx:) i — l XQ->--|-G<; 



where A'J^ stands for the axial- vector current. In the xQSM, the correlation function /C can be expressed as a functional 
integral 



/C ^ i / DijDij^DUJN{x)QHz)d^AliOJl{y)e^p 



cV'^ (ip + i^M{-d'^)W'- v/Af(-o»)2) V 



(10) 



in the chiral limit, where M(— 9^) denotes the momentum-dependent dynamical quark mass and U^'-" represents the 
chiral field defined as 

W^ = i±Ii[/ + l^^[/t (11) 

2 2 ^ ' 

with the Goldstone boson field U — exp(iA"7r''//^). 

It is, however, extremely complicated to solve Eq. (|10l) numerically, since the PV ttNN coupling constant involves the 
two-body quark operators Q* and the axial-vector one, which will lead to laborious triple sums in quark levels already 
at the leading order. Moreover, the momentum-dependent quark mass, which is known to be of great significance 
in describing nonleptonic processes ;27|, introduces in addition technical difficulties [30|. One way to avoid these 
complexities is to use a gradient expansion taking {pU /M) <^ 1 [3l| or equivalently is to start from the effective weak 
chiral Lagrangian already derived in Ref. j29]. Note that though we did not carry out the derivative expansion to 



order p'*, it is not difficult to estimate how large the corresponding LECs could be. In Ref. [28|, the AS' = 1 effective 
weak chiral Lagrangian to order p was investigated in the case of the local chiral quark model. As one can see, all 
of the LECs are order-of-magnitude smaller than the 0{p^) LECs. In this sense, even though we go further beyond 
the leading order, the contribution from higher derivative terms will not enhance or suppress h-j^NN rnuch. It will 
be at most below (5 — 10) %. Thus, we will use the AS" = effective weak chiral Lagrangian derived in Ref. f29l| as 



our starting point, instead of dealing with Eq. (|10l) . Nevertheless, the present approach goes beyond the previous 
analyses in the Skyrme model |19l - [2l| , because the present scheme incorporates properly the effects of the perturbative 
quark-gluon strong interaction in the derivation of the PV ttNN coupling constant. 

The leading-order (LO) term of the AS* = effective weak chiral Lagrangian in the large Nc can be expressed in 
terms of the vector and axial- vector currents 

2 

-Clo = 2( (5ii ^^ 

1=4 






(12) 



where the vector and axial- vector currents are defined as 

ya ^ IIt.it^^R^ + 7,^)] , ^a ^ IIt,[t^^R^ _ L,)] (13) 

in terms of L^ — iWd^U, R^ — iUd^U\ and T" — ( i ^, • • • , ^ ) . The parameter /^ stands for the pion decay 
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constant /^ — 93 MeV. The explicit expressions for the coefficients aij (a — a, (3, 7, p) can be found in Ref. [29^ . 
The classical soliton field Uq is assumed to have a structure of the trivial embedding of the SU(2) hedgehog field as 

^^^/exp(.r^fP(r)) 0\ ^^^^ 

with the profile function of the soliton P{r). This classical soliton field can be fiuctuated in such a way that the pion 
field can be coupled to a weak two-body operator 

U = exp(iT •7r/2)C/o exp(i-r-7r/2). (15) 

Similarly, the vector and the axial- vector currents transform as 

J n J TT 



A^^ a; + f /"''v>\ f; = f; + -r'^Ay, (le) 



where the indices a, b — I, ■ ■ ■ ,8 and i = 1, 2, 3. The current with a tilde indicates that arising from the background 
soliton field. 

Since the PV ttNN interaction Lagrangian is expressed as 

-Cpv = --=/^ijvw*jv(T X 7r)3*jv , (17) 

which is linear in the pion field and defined in the SU(2) flavor space (proportional to (r x 7r)3), one can easily see 
that the term X]i=4 ^^^ does not contribute to the PV nNN Lagrangian. Moreover, since /^^' — /°''* = 0, the pion 
fields for the PV ttNN Lagrangian can survive in the vector and axial- vector currents only when a = i. Writing them 
explicitly, we have 

^; - 4 + ^(^M X ^)^ ^' = V; + Ua, X ny. (18) 

Jtt Jtt 

Considering the terms contributing to the PV wNN vertex, we obtain for the LO Lagrangian 



ao = 511 E K A'^ + (3712 - 722) V^ AO^ + (3721 - 722) V^ A^^^ 



i=l 

+2722 V^ A^'' + ^722 {V^ A'"' + V^ A^^) + {V^ A). (19) 

v3 

Extracting the terms linear in the pion field from Eq. ([T^ and rearranging them, we finally derive the LO PV ttNN 
Lagrangian; 

In 

+ (3721 - 722)"^;°(V^'' X 7r)=' + (3712 - 722)A° (A^ X nf 

+-^722 [AliA^ X Tvf + V^iV^" X 7r)3] | + ('o"'^'^ O (O x tv)A . (20) 



For simplicity, we have omitted the tildes in the currents. 

In a similar manner, the next-to-leading order (NLO) effective weak chiral Lagrangian in the large N,. expansion 
derived in [23| yields the Lagrangian for the PV t:NN vertex as 
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(an -I- 2/3ii)(A3 x 77)3 + (^22 + 2/322) 



(A4 X 7r)4 + (A5 X 7V) 



4X1J3 

-^2 



+ 1 (7i2 + 2/5i2)(AoX7r)3 + 3 



I ^1 



- 1 (721 +2/52i)(Ao X 7r)3 



-2(722 -I- 2^22) 



(A3 X tt)3 + ^{As X 7r)3 



where 



J 7T 



Tr[{R^XaR^ + L^XaL'')T] for a = 0, 3, 



(A, X Tv)a = ^Tr [{R^XaR'' + L^XaL^')faMXbTr^] for a = 4, 5, 



(21) 

(22) 
(23) 



and Aq is defined as the unit matrix in SU(3) divided by 3. The integrals Xi in Eq. (j2ip were already evaluated in 
Ref. [29i] and are expressed as 



2:1 = - 
X2 = 
^3 = 



d^k 



M{k) 



(^^) 



M 



(27r)4 fc2-|-M2(fc) ANc ' 

d^fc M^{k)~!^M{k)M' _ fl 

(27r)4 (fc2 + Af2(fc))2 ~ 4iVe 



df^k 



^M"k'^ + \M' - |^fc2 M + M'^M' + ^M'^M" -t- ^k^MM'"^ + ^M' 



r 1 A>"1,2 1 1 A>' _ M!± 1.2 



fc2 -I- M2(fc) 

iM -f 2M2Af' -I- M3M'2 



(fc2+Af2(fc))2 



(fc2+Af2(fc))3 



where {ipip)M denotes the quark condensate in Minkowski space and M' = {dM{k)/dk)/2k. 
The next step is to carry out the zero-mode collective quantization of the soliton 

Ua{x) ^ U{x,t) ^ R{t)Uo{x)R\t) , 



(24) 
(25) 

(26) 

(27) 



where R{t) stands for the unitary time-dependent SU(3) orientation matrix of the soliton R{t) = exp{iil°-{t)X"' /2) 
with its angular velocity il°-{t) that is of order 0{1/Nc). Each current is transformed as 



J TT r 



A-^ 



Vo' = f^r{^-R[[Uo,R^R],U^,]R^ 



^ [1 - cos P{r)]D''°'n°' + fl sin^ P{r)D''^Ve - fl sin^ P{r){r ■ f2)£)" 



J_K_r 

2 



Vf = ^Tr{z'-R[Uo,m^o]R^]=ff-^^^^^^n,kfjD'^' 






^a ^ _:|xr('zyi?{[C/o,i?U],(7(l}i?t 



/2 f sin P{r) cos P{r)f'e,,kD'"'n^ + sm P{ry f^^^D'^'^nA 

in 



A- = :^Tr(*^i?{C/o,a.f/o^}i?n=/2 



sin2P(r) ^ 
2r '^ 



,^2_sin2PW 
2r 



Ijaj 



(28) 
(29) 

(30) 

(31) 



where Italic (Greek) indices run over 1, 2, 3 (4, • • • , 7), respectively, and dot (prime) means the derivative with respect 
to time (radius), respectively. The Wigner D functions and the angular velocity are defined as 



D''''{R) = iTr(A"i?A''i?^), R^ R = -A"f^°. 



(32) 



Before we proceed the calculation of h\j^^, we want to emphasize that we will investigate h\j^j^ first in the SU(2) 
case in this work. Of course, the strange quarks may still play a certain role in describing h\j^j~^- In fact, Ref. j35l | 
showed that the strange quark operator {q}?i,iQ){s"iin?<s ) in duced by Z° exchange could contribute significantly to 
the NN coupling constant. The main argument of Ref. |35l| lies in the fact that the A/ — 1 operator proportional 
to h^pfj^ can be related to the AS — 1 operator by an SU(3) rotation followed by an isospin rotation. Then, it was 
found that the linear combination of the strange operators made a large contribution to h^j^pj , which indicates that 
it has large matrix elements in the nucleon state. The SU(3) Skyrme model came to the similar conclusion that the 
four quark operators with the strange quark contributed to /i^jya' significantly [21| because of the induced kaon field. 
Note, however, that Ref. [21| has not used the renormalized effective weak Hamiltonian but started from the bare 
Hamiltonian. On the other hand, in a recent lattice study [38], the strange quark operators can only contribute to 
the quark-loop diagrams for which the signal-to-noise ratio remains far too small to bring out any reasonable signal, 
so that they were neglected. Moreover, recent findings have it that the content of strange quarks in the nucleon in the 
vector channel is neg ligible small 36| and that the strangeness in the scalar and axial- vector channels is still hampered 
by uncertainties [37|. Thus, it is still too early to reach a conclusion on the contribution of strange quark operators 
to ft-^Tvjv- Ii^ *he present work, we will concentrate on the case of SU(2), since it does not vanish even in SU(2). As 
we will discuss later in detail, this finite result is distinguished from that of the SU(2) Skyrme model [34] in which 
turns out to be equal to zero. The extension of the investigation to SU(3) will be found elsewhere, 
e we will calculate /i^jviv i^ ^he process of n7r+ — >■ p, we can rewrite the LO and NLO Lagrangians in SU(2) as 
follows: 

^Lo = ^^ {(-«ii + 2722) [V^V+^' + AlA+f^] + 2(3721 - 722)V;°y+'^ 
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/ 411X3 
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- 1 (721 + 



2/52i)A^ 



-^2 



2(72 



1 ) (712 + 2/5i2)A+ 
2p22)(A+ + A+) 



(33) 



(34) 



where we have used the identity 

(O X 7r)3 = V2i(7r"0+ - tt+O") 

with the definitions O^ = \{0'^ ± iO^) and vr* 

becomes the unity matrix with factor I/a/S in going from SU(3) to SU(2). The PV nNN coupling constant, h]^j^j^, 
can be directly read from the matrix element 



(35) 
4=(7ri ± iTT2). The eighth component of the Gell-Mann matrices 



hi 



NN 



i{pt\C^v\nt,7^+) , 



(36) 



where 






'^■nNN" 



f*jvi(7r r+-7r+r )*jv- (37) 

Let us first compute h^j^j^ with the LO Lagrangian. Note that the iso-scalar current vanishes identically in the 
present model. By using the results in the previous Section, we can see that the temporal component can contribute 
to /i^rATTv because of the orthogonality of D"'^. This produces the following expression: 



v,X^ = 






iD^^jn^ + D^'n^iD'^' + iD'^')n^ + D^'Vt^{D^^ + iD^^)Vl' 



"^'^^ - 3^/2 



sin^ P{r) cos^ P{r) 



D^'n^{D^' + iD^')n^ - D^^n^{D^3 + iD^^)n' 



(38) 



Because of the zero- mode quantization, the angular velocity is expressed in terms of the spin operator S'* 17* = S^ /I, 
where / is the moment of inertia of the soliton. The spin operator and Wigner D function satisfy the commutation 
relation [S\ D°-^] = ie^^^D"^'^ . Then, the matrix elements of Ea. (l38|) are written as 



{p t \V,X^\n t) = J^ ^ sin^ P{r) = -{pt \V+Vi 



{pt\AlA+\nt) 



15%/2 2/2 

JtA_ 

3^2 2/2 



l^t), 



sin^ P{r) cos^ P{r) = - (p ^ \A+ Al\n t) ■ 



(39) 



Since the LO Lagrangian is symmetric under the exchange of the indices 3 and +, it turns out that 

hlj,^iLO) = . (40) 



This null result of the LO h\j^j^ was also obtained in the minimal Skyrme model 34 1 



The NLO Lagrangian has a rather complicated structure, so that it is convenient to analyze first Aq 3. Introducing 



r], and I], as 



R, = -r%, L, = -T%, (41) 

we rewrite the expressions for A*- as 



K = YTr[(i?^T,i?^ + L^TaL'^y] = ^{r'^r"^ + ;™Z"^)Tr(T"V„r"T^) 
= y Cra^r^'' - Sy^r^'' + r'^Ta^, + {r ^ l)\ , for a ^ , 

K = Y|Tr[(i?^i?^ + L^L^)t'] = Y|(r;"r"'^ + /™r'^)Tr(r'"T"r') 



(42) 



Here, index a runs over a = 1, 2, 3. Then, r* and ^* become 

rg = -i:)'^''(-sinF(r)cosF(r)e^'™f7'F" + sin2p(r)5|if70 , (44) 

bf^i^bi. D/ N , cb,sin2P(r) sin^ F(r) ^fc, ^ 



Y 2r r y 

/g = -D"^(sinP(r)cosP(r)e'''"17'f" + sin2p(r)4'0') , (46) 

/," = -i^'^^f - f^r'drPir) - 4^!i^i^ + 5i!!!ZM^...,A (47) 

\ 2r r J 

where the transverse Kronecker delta is expressed as 5?^ — 6"^ — f°-f^ . Putting these results together, we arrive at 
the expressions for A'^ and Ajj": 



A+ = :|Tr[(P^r3P^ + L,T:,L^)r+] = ^ 



Url + iriy^ + rl^'' + i^) + (r ^ /) j , 



A+ = :|Tr[(P^AoP'^ + L^AoP^)r+] = ^ (^(r^ + triy^ - riyf^ + i^) + (r ^ /)) . (48) 

Since 

f d^x{p-t\rl{r^''+ir^^)\nt) = - f d^x {p t \ir^^ + ir^^)rl\n t) 

= -^ I drr^ sin^ P{r) {sin^ P(r)-3cos^P{r)), (49) 

one can easily see that only A^ contributes to h\^j^. As a result, h\j^j^ from the NLO Lagrangian turns out to be 
/lijvw(NLO) - 1^ ('AA9 + ^AAio')ydrr2sin2p(r)(sin2p(r)-3cos2p(r)), (50) 

where the LECs A/g and A/io are given as |23l 

AAg = 4iVc [4Iil3(7l2 + +721 + 2/5i2 + 2p2l) + I|(712 - 721 + 2pi2 - 2p2l)] 

/4 

= 4(-0'0)Af2:3 (712 + 721 + 2pi2 + 2p2i) + j^ (712 - 721 + 2(012 - 2p2i) , (51) 

Nw = 4iV,Z| (722 + 2^22) = |^(722+/522). (52) 

As we will discuss later, the LECs Ag and A/10 are essential to describe the PV ttNN coupling constant. 



III. RESULTS AND DISCUSSION 



We are now in a position to calculate Ea.(|50p numerically. In doing so, we make use of the momcntumc dependent 
quark mass derived from the instanton vacuum [26| and the corresponding results of the LECs obtained in Ref. [2^ . 
The value of Mq = M{k = 0) is taken to be 350 MeV as in Ref. [29|, which was fixed by the saddle-point equation 
from the instanton vacuum j26| . Moreover, we employ three different types of the solitonic profile function to examine 
the dependence of /i^^jy '^"^ them. The first one is the arctangent profile function P{r) [32l | 



P(r) ==2arctan('— ] , 



(53) 



where rg is given by tq 






Employing qa — 1.26 and /^r = 93 MeV, we obtain tq = 0.582 fm. We use a 



physical profile function as a second one, which associates with the proper pion tail of the nucleon 

pi^\ = / 2 arctan (^) (r < r^^) 

^' \ Ae-"''''(l + m^r)/r2 {r > r.^), 



(54) 



where ttIt^ denotes the pion mass and A — 2r§. r^, is determined by the intersection of the arctangent function (r < Vx) 
and pion tail (r > Tx). If one takes the limit ttItt — >■ for the pion tail, the physical profile function becomes identical 
with the arctangent one at large r. With the physical pion mass considered, we have r^ = 0.749 fm. The final one is 
the linear profile function initially proposed by Skyrme [33| 






./A) 



(m< A) 
{u>X) 



(55) 



where u = 2ef-^r with e — 4.84 and A = 3.342. 

Using these three profile functions, we can immediately compute the PV ttNN coupling constant h\^j^. Figure [1] 
draws the results of /i^jvat a-s a function of the Wilson coefficient K. The solid curve depicts that with the physical 




K 



FIG. 1: (Color online) PV ttNN coupling constant /i^jvjv as a function of the Wilson coefficient K in units of 10~*. The 
solid curve draws the result with the physical profile function, the dashed one depicts that with the arctangent profile, and the 
short-dashed one does that with the linear one. 



profile function, whereas the dashed and short-dashed ones correspond to those with the arctangent and linear profile 
functions respectively. One can regard the difference between the results with the physical profile function and those 



with the arctangent one as effects of the finite pion mass, which contribute to h^j^^ approximately by 10 %. Moreover, 
the type of the profile function does not change much the general features of h^f^j^, though we preferably take the 
results with the physical one as our final values. 

We find out from Fig. [T] that h\j^j^ is rather sensitive to the Wilson coefficient K and it decreases monotonically, 
as K increases. Wc notice that its sign is even changed around K — Q. This can be easily understood. The LECs 
A/g and A/io in Eq. (|50|) play essential roles in determining the K dependence of h\jqj^. Figure [2] draws the results of 
the LECs A/g and A/io as functions of K. While Ag depends rather strongly on K , JVio does mildly on K. Moreover, 
A/g is dominant over A/io, so that the PV nNN coupling constant is mainly governed by A/g. Since Ag is the main 
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FIG. 2: (Color online) Low-energy constants Ag, A/io, and J\fg + 2A/io/3 as functions of the Wilson coefficient K in units of 
lO"'^^ GeV'^. The solid curve draws the result of Ag + 2A/io/3, dashed one depicts that of A/g, and the short-dashed one does 
that of A/io- 



contribution to h\j^pf, we want to examine it in detail. We can easily see that the first term of Ea. (j5ip containing the 
quark condensate is much larger than the second one. Moreover, since 712 and p2i are much smaller than the other 
two coefficients 721 and pi2, we can neglect them in Ag. Then, A/g can be expressed as 



A/'g « 4(i/'V')m2'3 (721 + 2pi2) 



(56) 



As shown in Eq.([T]), the coefficient 721 comes from the original effective weak Hamiltonian at the mass scale of the W 
boson fx = Mw — 80.4 GeV corresponding io K = 1. In this case, only 721 survives in A/g. However, when we start to 
scale the Hamiltonian down to /i ~ 1 GeV that corresponds to -ftT ~ 4, the gluonic renormalization arising from giuon 
exchange parallel to Z-boson exchange is turned on. As a result, the 2pi2 term becomes as large as a half of the 721 
one [3, [29| at this scale. If one goes further down to the scale at which iiT « 6 ^, the correction of pi2 cancels out the 
contribution of 721, so that h^pfpf almost vanishes, as already shown in Figs. [1] and [2j This cancellation implies that 
the effects of the gluon renormalization leads to the suppression of the PV ttNN coupling constant in the present 
approach of the SU(2) xQSM. 

Since Ref. [23| derived the effective weak chiral Lagrangian based on the AS = effective weak Hamiltonian, it is 
plausible to take the value K = A for h^pfM, which corresponds to the renormalization scale of the charm quark mass 
II ~ I GeV as done for the AS = 1 case Q- Thus, we obtain h]^pfj^ « 1 x 10^* for K = A. However, if one neglects 



^ There is a caveat in scaling further down below 1 GeV, because the matching problem becomes non-trivial below the charm quark mass. 
Furthermore, we still do not know how to incorporate all possible nonperturbative effects consistently below 1 GeV. 
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all the renormalization effects, i.e. if one takes K ^ 1, we have hl^jvM « 4 x 10~*, which is similar to that of the 
SU(2) Skyrme model with vector mesons {h^j^j^ = (2 — 3) x 10~^) [1^ in which the effective weak Hamiltonian at 
/i = Mw or with K = 1 was used. Note that the present result is almost 40 times smaller than the "best" value of 
DDH (/^i^j^ = 4.5 X 10-^) i. 
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FIG. 3: (Color online) PV nNN coupling constant as a function of the pion mass m,i 
iC" = 4 is used, which corresponds to /i ~ 1 GeV. 



in units of 10 . The Wilson coefficient 



Very recently, Ref. ^ has reported the first resuh of lattice QCD: h^j^j^ = (1.099 ± 0.505loo64) x 10"^ with the 
pion mass m^ = 389 MeV. Thus, it is interesting to compare the present result with the lattice one. In order to do 
that, we need to examine the dependence of /i^jvAf '^^ ^^^ pion mass TOtt. Figure |3] depicts the PV ttNN coupling 
constant as a function of to^. We employ here the physical profile function with K = 4. Interestingly, h\^j^ starts 
to increase, as m,^ does. As a result, ft-^^^v turns out to be around 1.8 x 10~® for m^ — 400 MeV. Though it is 
still around five times less than that of the lattice calculation, we can infer from Fig. [3] that lattice results with the 
physical pion mass might be quite smaller than that of Ref. 1381 ]. Moreover, if one takes K — 1, /i^jvat with ra-^ = 389 



MeV would become hi 



■NN 



6.77 X 10 that is comparable to the lattice one, though the value K = 1 does not seem 



tenable for h^j^j^ as discussed previously. However, one has to keep in mind that Ref. [38] has not performed the 
calculation of the full matrix element, since the quark-loop diagrams were omitted because of technical difficulties. A 
quantitative comparison with full lattice calculations is still being awaited. 



IV. SUMMARY AND OUTLOOK 



We have investigated the parity- violating pion-nucleon coupling constant /i^jvat within the framework of the chiral 
quark-soliton model with the gradient expansion used. Starting from the AS* = effective weak chiral Lagrangian 
derived in the same framework, we have calculated the parity- violating ttNN coupling constant. It was found that it 
vanished at the leading order in the large Nc, i.e. h]^j^j^(LO) = 0, which is of order 0{Nc 



-l/2\ 



but it was finite to 



the next-to-leading order, i.e. 0{Nc )■ 

Employing three different profile functions, that is, the arctangent, physical, and linear ones, we calculated the 
parity- violating -nNN coupling constant to the next to the leading oder. It turns out that the values of h\j^p^ depend 
sensitively on the values of the Wilson coefficient K and vanishes around K rv Q. The reason can be found in the fact 
that the contribution of the gluonic renormalization constant /5i2 cancels out that of the 721, which is the leading one. 
It indicates that the perturbative gluonic contribution supresses the parity- violating ttNN coupling constant. Taking 
the scale of the charm quark mass, i.e. /x « 1 GeV, we found h^j^j^ « 1 x 10~^, which is almost 40 times smaller 
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than the "best vahie" of Ref. Q . If the u = Mw is selected, the value of h\j^^ turns out to be similar to that from 
the Skyrme model with vector mesons [1^ . We also compared the present result with that of the lattice calculation. 
Thus, we examined the dependence of the parity-violating TriViV coupling constant on the pion mass and found that 
^nNN increased as to^ did. If one uses tti^ = 400 MeV, the result turns out to be almost two times larger than that 
with the physical value m,^ = 140 MeV but is still about five times smaller than the lattice one. However, we want to 
emphasize that neither the present result nor the lattice one is the final one. 

In order to understand the parity- violating ttNN coupling constant h^j^pj more completely and quantitatively, we 
have to consider the following important physics: Since we have considered the SU(2) case in the present work, the 
effects of strangeness were left out. As already mentioned in Section |TT] , however, the strange quark operators may 
play a certain role in describing the parity-violating ttNN coupling constant. Extending from SU(2) to SU(3) is 
lengthy but straightforward in the present framework. Starting from Eqs. (|20| I2ip . we employ the quantization with 
the embedding p^ . In particular, the singlet current is distinguishable from the octet one in SU(3), so that this 
would make difference in predicting the parity-violating ttNN coupling constant. Moreover, the fourth and fifth flavor 
components of the vector currents enter the next-to-leading order Lagrangian, this would also contribute to /i^^vjv ■ 

As was seen in Section [nil the gluon renormalization plays a role of suppressing the parity-violating irNN coupling 
constant. However, the effective weak Hamiltonian at two- loop order was derived very recently in Ref. [39|, where 
the QCD penguin diagrams were also considered. This Hamiltonian is more complete than that from Ref. [8|. Thus, 
it is of great significance to investigate the AS* = effective weak chiral Lagrangian based on this effective weak 
Hamiltonian. It is also of great interest to see whether these penguin diagrams enhance the h\jqj^ or not, since they 
give part of answers of explaining the A/ = 1/2 rule in nonleptonic decays [y, 0]- The corresponding investigations 
are under way. 
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